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ABSTRACT 

The  second  part  of  this  work  investigates  the  seasonal  variabilities  of  the  Japan/East  Sea  (JES)  circulation 
using  the  U.S.  Navy  Generalized  Digital  Environmental  Model  (GDEM)  climatological  temperature  and  salinity 
dataset  (public  domain)  on  a  0.5°  X  0.5°  grid.  A  variational  P-vector  method  was  developed  to  invert  the  velocity 
field.  The  GDEM  for  the  JES  was  built  up  on  historical  (1930-97)  136  509  temperature  and  52  572  salinity 
profiles.  The  climatological  mean  and  seasonal  variability  of  the  current  systems  are  well  inverted,  especially 
the  Tsushima  Warm  Current  and  its  bifurcation,  the  East  Korean  Warm  Current  (EKWC),  the  Japan  nearshore 
branch,  the  confluence  of  the  EKWC,  and  the  North  Korean  Cold  Current  near  the  Korean  coast  and  flows 
northeastward  along  the  subpolar  front,  and  a  mesoscale  anticyclonic  eddy  in  the  Ulleng/Tsushima  Basin.  Fur¬ 
thermore,  this  method  has  the  capability  to  invert  flow  reasonably  well  across  the  shallow  straits  such  as  the 
Tsushima/Korea,  Tsugaru,  and  Soya  Straits.  The  GDEM  temperature  and  salinity  and  the  inverted  velocity  fields 
provide  balanced  initial  fields  for  JES  numerical  modeling  and  simulation. 


1.  Introduction 

The  Japan  Sea,  known  as  the  East  Sea  in  Korea,  has 
a  steep  bottom  topography  (Fig.  1)  that  makes  it  a  unique 
semi-enclosed  ocean  basin  overlaid  by  a  pronounced 
monsoon  surface  wind.  The  Japan/East  Sea  (JES)  covers 
an  area  of  106  km2,  has  a  maximum  depth  in  excess  of 
3700  m,  and  is  isolated  from  open  oceans  except  for 
small  (narrow  and  shallow)  straits  that  connect  the  JES 
to  the  Pacific  Ocean.  The  JES  contains  three  major  ba¬ 
sins  called  the  Japan  Basin  (JB),  Ulleng/Tsushima  Basin 
(UTB),  and  Yamato  Basin  (YB),  and  a  high  central  sea¬ 
mount  called  the  Yamato  Rise  (YR).  The  JES  has  a  great 
scientific  interest  as  a  miniature  prototype  ocean.  Its 
basinwide  circulation  pattern,  boundary  currents.  Sub¬ 
polar  Front  (SPF),  mesoscale  eddy  activities  and  deep 
water  formation  are  similar  to  those  in  a  large  ocean. 

The  JES  experiences  two  monsoons,  winter  and  sum¬ 
mer,  every  year.  During  the  winter  monsoon  season,  a 
cold  northwest  wind  blows  over  the  JES  as  a  result  of 
the  Siberian  High  Pressure  System  located  over  the  East 
Asian  continent.  Radiative  cooling  and  persistent  cold 
air  advection  maintain  cold  air  over  the  JES.  The  north- 
west-southeast  oriented  jet  stream  is  positioned  at  the 
JES.  A  typical  winter  monsoon  pattern  lasts  nearly  six 
months  (Nov  to  Apr).  During  the  summer  monsoon,  a 
warm  and  weaker  southeast  wind  blows  over  the  JES. 
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A  typical  summer  monsoon  pattern  lasts  nearly  four 
months  (mid-May  to  mid-Sep). 

The  JES  general  circulation  has  been  investigated  for 
several  decades.  The  Tsushima  Warm  Current  (TWC), 
dominating  the  surface  layer,  flows  in  from  the  Tsushima 
Strait,  and  carries  warm  water  from  the  south  up  to  40°N 
where  a  polar  front  forms  (Seung  and  Yoon  1995).  Most 
of  the  nearly  homogeneous  water  in  the  deep  part  of 
the  basin  is  called  the  Japan  Sea  Proper  Water  (Moriyasu 
1972)  and  is  of  low  temperature  and  low  salinity.  Above 
the  proper  water,  warm  and  saline  water  that  enters 
through  the  Tsushima  Strait  flows  northeastward  and 
flows  out  through  the  Tsugaru  and  Soya  Straits. 

The  TWC  separates  north  of  35°N  into  western  and 
eastern  channels  (Uda  1934;  Kawabe  1982a, b;  Chu  et 
al.  2001a;  Chu  et  al.  2001b).  The  flow  through  the  west¬ 
ern  channel  closely  follows  the  Korean  coast  [called  the 
East  Korean  Warm  Current  (EKWC)]  until  it  bifurcates 
into  two  branches  near  37°N.  The  eastern  branch  follows 
the  SPF  to  the  western  coast  of  Hokkaido  Island,  and 
the  western  branch  moves  northward  and  forms  a  cy¬ 
clonic  eddy  at  the  Eastern  Korean  Bay  (EKB).  The  flow 
through  the  eastern  channel  follows  the  Japanese  coast, 
called  the  “Nearshore  Branch”  by  Yoon  (1982).  More 
accurately,  we  may  call  it  the  Japan  nearshore  branch 
(JNB).  The  JNB  is  usually  weaker  than  the  EKWC.  The 
TWC  at  both  channels  reduces  with  depth.  The  North 
Korea  Cold  Current  (NKCC)  meets  the  EKWC  at  about 
38°N  with  some  seasonal  meridional  migration.  After 
separation  from  the  coast,  the  NKCC  and  the  EKWC 
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Fig.  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the  Japan/East  Sea 

(JES). 


converge  and  form  a  strong  front  (i.e.,  SPF)  that  stretch¬ 
es  in  a  west-east  direction  across  the  basin.  The  NKCC 
makes  a  cyclonic  recirculation  gyre  in  the  north,  but 
most  of  the  EKWC  flows  out  through  the  outlets  (Uda 
1934).  The  formation  of  NKCC  and  the  separation  of 
EKWC  are  due  to  a  local  forcing  by  wind  and  buoyancy 
flux  (Seung  1992).  Large  meanders  develop  along  the 
front  and  are  associated  with  warm  and  cold  eddies. 

The  seasonal  variability  of  the  TWC  at  the  Tsushima/ 
Korea  Strait  largely  impacts  the  JES  physical  condi¬ 
tions.  An  accurate  estimate  of  the  volume  transport  at 
this  strait  is  important  for  the  JES  circulation  and  ther¬ 
mohaline  structure.  Two  methods  have  been  used  to 
estimate  the  volume  transport  through  the  Tsushima/ 
Korea  Strait:  1)  dynamical  calculation  of  the  hydro- 
graphic  data  with  a  level  of  “no  motion”  at  the  bottom 
and  2)  calculation  of  the  sea  level  difference  across  the 
strait. 

Using  the  dynamical  calculation  approach,  the  vol¬ 
ume  transport  through  the  Tsushima/Korea  Strait  is 
characterized  by  the  minimum  in  winter-spring,  and 
maximum  in  summer-fall  (e.g.,  Hidaka  and  Suzuki 
1950;  Yi  1966).  For  example,  Yi  (1966)  estimated  the 


annual  average  of  1.3  Sv,  a  maximum  of  2.2  Sv  in 
October,  and  a  minimum  of  0.3  Sv  in  February  (Sv  = 
106  m3  s~')-  Later  on,  Mitta  and  Ogawa  (1984)  analyzed 
the  historical  velocity  data  obtained  by  the  current  mea¬ 
surements  lasting  one  day  in  the  summers  of  1942  and 
1943  to  investigate  the  current  structure  across  the 
Tsushima/Korea  Strait.  The  data  show  intense  north¬ 
ward  currents  near  the  bottom  in  the  Tsushima/Korea 
Strait.  This  leads  to  doubt  about  the  validity  of  the  dy¬ 
namical  calculation  with  the  level  of  no  motion  at  the 
bottom  of  the  strait  (Mitta  and  Ogawa  1984).  To  remove 
the  “level  of  no  motion”  assumption,  we  use  an  inverse 
method  to  obtain  absolute  geostrophic  velocity  from 
hydrographic  data. 

Using  the  sea  level  difference  across  the  strait  without 
considering  the  contribution  from  baroclinic  motion, 
Kawabe  (1982a)  and  Toba  et  al.  (1982)  found  the  vol¬ 
ume  transport  through  the  Tsushima/Korea  Strait  to  be 
the  same  as  the  dynamical  calculation  (Hidaka  and  Su¬ 
zuki  1950;  Yi  1966).  Considering  the  effect  of  baroclinic 
motion  and  subtracting  the  sea  level  difference  due  to 
the  baroclinic  motion  from  the  observed  current  data 
(1988-90),  Isobe  (1994)  found  that  the  volume  transport 
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Table  1.  Glossary  of  acronyms. 


Acronym 

Explanation 

EKB 

Eastern  Korean  Bay 

EKWC 

East  Korean  Warm  Current 

GDEM 

U.S.  Navy  Generalized  Digital  Environmental  Model 

JB 

Japan  Basin 

JES 

Japan/East  Sea 

JNB 

Japan  nearshore  branch 

MOODS 

U.S.  Navy  Master  Oceanographic  Observational 
Data  Set 

NKCC 

North  Korean  Cold  Current 

SPF 

Subpolar  Front 

TWC 

Tsushima  Warm  Current 

UTB 

Ulleng/Tsushima  Basin 

YB 

Yamato  Basin 

YR 

Yamato  Rise 

was  at  a  maximum  in  early  winter  and  a  minimum  in 
early  spring  with  the  annual  range  (maximum  minus 
minimum)  of  the  seasonal  variation  as  0.7  Sv. 

Most  of  these  analyses  of  seasonal  variability  of  the 
JES  circulation  and  the  volume  transport  at  the  Tsush¬ 
ima/Korea  Strait  were  based  on  temporally  and  spatially 
limited  data.  Using  a  more  complete  dataset,  we  may 
improve  these  analyses  and  get  statistically  significant 
results.  The  U.S.  Navy  Master  Oceanographic  Obser¬ 
vational  Data  Set  (MOODS)  contains  136  509  temper¬ 
ature  and  52  572  salinity  profiles  during  1930-97.  Based 
on  the  MOODS  data,  an  unclassified  Generalized  Digital 
Environmental  Model  (GDEM)  was  established  with 
climatological  annual  and  monthly  mean  temperature 
and  salinity  fields  on  a  0.5°  X  0.5°  grid. 

Chu  et  al.  (2001a)  used  the  P-vector  inverse  method 
(Chu  1995,  2000;  Chu  et  al.  1998a,b)  to  calculate  the 
JES  absolute  geostrophic  velocity  from  the  GDEM  data. 
The  P-vector  method  contains  two  steps:  1)  determi¬ 
nation  of  the  velocity  direction  and  2)  determination  of 
the  velocity  magnitude.  Two  necessary  conditions  for 
the  inversion  are  easily  implemented  into  this  method: 
1)  the  isopycnal  surface  does  not  parallel  the  isosurface 
of  potential  vorticity  and  2)  the  velocity  has  a  vertical 
spiral  (Chu  et  al.  1998a,b).  The  P-vector  method  inverts 
the  circulation  in  the  JES  basin  reasonably  well  (Chu 
et  al.  2001a),  but  it  fails  to  invert  the  velocity  field  at 
shallow  straits  such  as  the  Tsushima/Korea,  Tsugaru, 
and  Soya  Straits. 

In  this  study,  we  develop  a  variational  P-vector  meth¬ 
od  to  improve  the  inversion  and  to  obtain  realistic  cir¬ 
culation  in  the  JES  basin  as  well  as  in  the  shallow  straits. 
The  outline  of  this  paper  is  as  follows.  A  description  of 
the  GDEM  data  is  given  in  section  2.  The  P-vector 
method  and  its  deficiencies  are  presented  in  section  3.  The 
variational  P  vector  is  presented  in  section  4.  The  inverted 
absolute  geostrophic  velocity  field  is  discussed  in  section 
5.  The  volume  transports  in  the  Tsushima/Korea,  Tsugaru, 
and  Soya  Straits  are  presented  in  section  6.  In  section  7 
we  present  our  conclusions.  There  are  quite  a  few  acro¬ 
nyms  used  in  this  paper,  listed  in  Table  1. 


2.  The  Navy’s  GDEM  Dataset 

Data  for  building  the  current  version  of  GDEM  cli¬ 
matology  for  the  JES  were  obtained  from  the  MOODS 
data.  MOODS  is  a  compilation  of  ocean  data  observed 
worldwide  consisting  of  1)  temperature-only  profiles,  2) 
both  temperature  and  salinity  profiles,  3)  sound-speed 
profiles,  and  4)  surface  temperature  (drifting  buoy).  The 
main  limitation  of  the  MOODS  data  is  its  irregular  dis¬ 
tribution  in  time  and  space.  Certain  periods  and  areas 
are  oversampled,  while  others  lack  enough  observations 
to  gain  any  meaningful  insights  (Chu  et  al.  1997a, b). 
Vertical  resolution  and  data  quality  are  also  highly  var¬ 
iable  depending  on  instrument  type  and  sampling  ex¬ 
pertise.  The  monthly  distributions  of  the  JES  tempera¬ 
ture  (Fig.  2a)  and  salinity  (Fig.  2b)  stations  show  that 
the  number  of  temperature  stations  is  2-3  times  more 
than  the  number  of  salinity  stations,  and  January  has 
the  least  profiles  and  August  the  most.  Yearly  temper¬ 
ature  (Fig.  3a)  and  salinity  (Fig.  3b)  profile  numbers 
show  temporally  uneven  distribution  with  almost  no  ob¬ 
servations  in  the  whole  JES  in  certain  years  (e.g.,  1944, 
1989  for  temperature,  and  1944,  1987-93  for  salinity) 
and  many  observations  in  other  years  (e.g.,  nearly  6500 
temperature  profiles  in  1969,  and  3700  salinity  profiles 
in  1967).  Spatial  and  temporal  irregularities  along  with 
lack  of  data  in  certain  regions  must  be  carefully  weight¬ 
ed  in  order  to  avoid  statistically  induced  variability. 
Based  on  the  MOODS  data,  the  Navy’s  global  clima¬ 
tological  monthly  mean  temperature  and  salinity  dataset 
was  built  (i.e.,  GDEM)  with  a  four-dimensional  (lati¬ 
tude,  longitude,  depth,  and  month)  display. 

The  basic  design  concept  of  GDEM  is  the  determi¬ 
nation  of  a  set  of  analytical  curves  that  represent  the 
mean  vertical  distributions  of  temperature  and  salinity 
for  grid  cells  (0.5°  X  0.5°)  through  the  averaging  of  the 
coefficients  for  the  curves  found  for  individual  profiles 
(Teague  et  al.  1990).  Different  families  of  representative 
curves  have  been  chosen  for  shallow,  middepth,  and 
deep-depth  ranges,  with  each  chosen  so  that  the  number 
of  parameters  required  to  yield  a  smooth,  mean  profile 
over  the  range  was  minimized.  As  pointed  out  by  Teague 
et  al.  (1990),  large-scale  oceanographic  features  are  gen¬ 
erally  found  to  be  similarly  represented  in  both  GDEM 
and  the  NOAA  Climatological  Atlas  of  the  World  Ocean 
temperature  and  salinity.  GDEM  appears  to  render  better 
representations  of  seasonal  variability  and  regions  of 
high  current  shear  because  of  a  different  smoothing 
method  and  a  finer  grid  spacing.  GDEM  data  contains 
the  monthly  mean  temperature  and  salinity  ( T ,  S )  and 
annual  mean  temperature  and  salinity  (T ,  S)  fields.  In¬ 
terested  readers  are  referred  to  Teague  et  al.  (1990)  for 
more  information. 

3.  P-vector  inverse  method  and  its  deficiency 

a.  Reduced  physics 

As  pointed  out  by  Wunsch  and  Grant  (1982),  in  de¬ 
termining  large-scale  circulation  from  hydrographic 
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data,  we  can  be  reasonably  confident  of  the  assumptions 
of  geostrophic  balance,  mass  conservation,  adiabatic, 
and  no  major  cross-isopycnal  mixing  (except  for  water 
masses  in  contact  with  the  atmosphere).  Under  these 
conditions,  the  density  of  each  fluid  element  would  be 
conserved,  which  mathematically  is  given  by 

V  •  Vp  =  0,  (1) 

where  p  is  the  potential  density  and  V  =  ( u ,  v,  w)  is 
the  geostrophic  velocity.  The  conservation  of  potential 
vorticity  equation  can  be  obtained  by  differentiating  (1) 
with  respect  to  z,  using  geostrophic  and  hydrostatic  bal¬ 
ances,  and  including  the  latitudinal  variation  of  the  Cor¬ 
iolis  parameter  to  give 

V  •  Vq  =  0,  (2) 

where  q  =  fdp/dz .  Equations  (1)  and  (2)  imply  that  the 
velocity  V  is  parallel  to  Vq  X  Vp. 

b.  Necessary  conditions 

Stommel  and  Schott  (1977)  pointed  out  that  the  three- 
dimensional  velocity  field  cannot  be  determined  from 
the  density  field  alone  when  the  q  and  p  surfaces  co¬ 
incide.  The  first  necessary  condition  for  the  validity  of 
this  inverse  method  is  as  follows. 

Condition  1:  The  p  surface  is  not  parallel  to  the  q  sur¬ 
face,  which  mathematically  requires 

Vp  X  Vq  ¥=  0. 

Stommel  and  Schott  (1977)  further  pointed  out  that 
the  three-dimensional  velocity  field  cannot  be  deter¬ 
mined  from  the  p  field  alone  if  the  horizontal  velocity 
does  not  rotate  with  depth  (/ 3  spiral).  The  existence  of 
the  (3  spiral  is  the  second  necessary  condition. 


V  =  r(x,  y,  z)P,  (4) 

where  r  is  the  proportionality.  Applying  the  thermal 
wind  relation  at  two  different  depths  zk  and  z,„,  a  set  of 
algebraic  equations  for  determining  the  parameter  r  is 
obtained : 

r®P®  -  =  A ukm 

r®P®  -  =  Avkm,  (5) 

which  are  two  linear  algebraic  equations  for  r,k)  and  r(m) 
[ru)  =  r(x,  y,  z,.)].  Here 

<4,“-A  v-)~  Lit  <6) 

where  p  is  the  in  situ  water  density,  and  pu  is  the  char¬ 
acteristic  value  of  the  density. 

The  existence  of  a  solution  of  (5)  implies  a  nonzero 
determinant  of  the  coefficient  matrix  of  (5),  which  is 
necessary  condition  2.  This  determinant  is  the  sine  of 
the  vertical  turning  angle  between  P;®  and  P}®  (Chu 
2000;  Chu  et  al.  1998a,b;  2001a). 

For  water  columns  satisfying  the  two  necessary  con¬ 
ditions,  we  solve  (6)  to  obtain  r(k>  for  the  level  z.k.  There 
are  N  —  1  sets  (m  =  1,2 ,  k  —  1,  k  +  1,  •  •  • ,  N)  of 
equations  (5)  for  calculating  r®.  Here  N  is  the  total 
number  of  vertical  levels  of  the  water  column.  The  N 
—  1  sets  of  equations  are  compatible  under  the  thermal 
wind  constraint  and  should  provide  the  same  solution. 
However,  because  of  errors  in  measurements  (instru¬ 
mentation  errors)  and  computations  (truncation  errors), 
the  parameters  r®  may  vary  with  m.  A  least  squares 
error  algorithm  is  used  to  minimize  the  error. 


Condition  2:  The  velocity  ( u ,  v )  should  execute  a  [3 
spiral,  which  mathematically  requires  that,  for  at  least 
two  depths,  z  =  Zk  and  z  =  zm,  with  horizontal  velocities 
[w®,  u®]  and  [«<"■>,  «<■>], 

K®  U® 

7^  0. 

M<m)  V(® 

If  we  cannot  find  levels  zk  and  zm  such  that  the  necessary 
condition  2  is  satisfied,  the  inverse  method  will  fail  to 
get  velocity  in  that  water  column. 


c.  P-vector  method 


Consider  the  unit  vector  P  (Chu  1995),  defined  by 


Vp  X  Vi q 
[Vp  X  Vpl 


(3) 


The  existence  of  this  unit  vector  implies  nonzero  de¬ 
nominator  of  (3),  which  is  the  necessary  condition  1. 
The  velocity,  V  =  (m,  v,  w),  parallels  the  unit  vector  P, 


d.  Deficiency  of  the  current  P-vector  method 

Let  (u"‘\  v(P))  be  the  absolute  velocity  determined  by 
the  P-vector  method,  and  (U(P),  V'p>)  be  their  vertical 
integrations, 


(U<p\  V(p>)  =  I  (u,p\  v(p>)  dz.  (7) 

J -H 

Due  to  the  local  determination  of  the  absolute  geostrophic 
velocity,  the  present  P-vector  method  does  not  always 
guarantee  mass  conservation  over  a  domain  cr,  that  is, 


ff 

dU(P) 

+ 

<3V®> 

J  Jcr 

dx 

dy 

dx  dy  ¥=  0 


(8) 


is  possible.  Such  a  deficiency  largely  affects  the  quality 
of  the  inversion.  For  example,  the  TWC  in  the  shallow 
Tsushima/Korea  Strait  was  not  well  inverted  (Chu  et  al. 
2001a).  Therefore,  we  develop  a  variational  algorithm 
taking  into  account  the  mass  conservation. 
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Fig.  2.  Spatial  distributions  of  MOODS  stations  during  1930-97:  (a)  temperature  and  (b)  salinity. 
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Fig.  2.  ( Continued ) 
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Year 


Fig.  3.  Temporal  distributions  of  MOODS  stations  during  1930-97:  (a)  temperature  and 

(b)  salinity. 


4.  A  variational  algorithm 

a.  Formulation 


Let  ( U ,  V)  be  the  vertically  integrated  velocity  op¬ 
timally  determined  by  minimizing  the  following  func¬ 
tional  (called  the  cost  function), 

J(U,V)  =  ^  J  J  [(U  -  L FP))2  +  ( V  -  V'<'-))2]  dx  dy  (9) 

with  the  constraint  (mass  conservation) 


dU(p)  dV <P) 

- +  - —  0. 

dx  dy 


(10) 


This  problem  becomes  an  unconstrained  optimization 
using 


L(U,  V,  A)  =  J(U,  V) 


+ 


ff  y 

dum 

dV(p> 

1  1  A 

+ 

J  Jo- 

dx 

dy 

dx  dy, 


where  A  is  the  Lagrangian  parameter. 


(11) 


b.  Horizontal  discretization 

Let  the  dependent  variable  f(x,  y)  be  defined  on  the 
interval  0  <  x  <  Lx,  0  <  L  .  Use  a  uniform  grid, 

0  =  Xj  <  x2  <  ■  •  •  <  xNi  =  Lx,  and  0  =  yj  <  y2  < 
•  •  •  <  yN  =  Lv  with  grid  spacing  of  (Ax,  Ay),  as 
shown  in  Fig.  4a.  Here, 

Ax  =  xi+l  -  x,  =  L  JNX,  Ay  =  yJ+l  -  y}  =  Ly/Ny. 

Let  p,  A  be  evaluated  at  the  grid  point  ( i ,  j),  and  the 
integrated  velocity  components  U,  V,  U"‘>,  V"']  be  eval¬ 
uated  at  the  staggered  points,  respectively  (Fig.  4b).  The 
functional  (11)  is  discretized  as 

J  Nx~  1  Ny—1 

i  =  -n  mj  -  Ujpy  +  (V„  -  V<py]AxAy 

^  i=  1  j=  1 

J  Nx-  1  Ny-l 

+  r  2  2  A, .(<4  +  uu  ,  -  Ut  UJ  -  U,  w  ,)Ay 

^  l=  1  j=  1 

J  Nx—  1  Ny-l 

+  r  2  2  A„d4  +  V,_hj  -  V,,  ,  -  V,  ,  )A.v. 

^  i=l  7=1 

(12) 
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Fig.  4.  Staggered  grid  used  for  the  computation:  (a)  grid  distribution  and  (b)  staggered  grids  for  ( u ,  v)  and 

standard  grids  for  (A,  p). 


c.  Combined  local-global  determination 


Minimization  of  L  becomes  a  combination  of  local 
determination  at  the  staggered  grid  (i,  j )  for  velocity  ((3- 
spiral  approach) 

uy  -a  uf,  v,  -a  vf  (i3) 

and  global  determination  (box  model) 


+  Uu-'  -  U‘-'j  ~  U  '  <  1 

+  ±<yv  +  y. -  y.j ,  -  y, ,, ;)  ->  o,  <  i4> 


which  keeps  the  mass  conservation  at  any  box  centered 
at  the  nonstaggered  grid  for  p  and  A  (Fig.  4b).  Thus, 
this  variational  P-vector  method  can  be  treated  as  a  f3- 
spiral  box  model. 
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Table  2.  Inverted  monthly  variation  of  volume  transport  (Sv)  at  three  major  straits.  The  positive/negative  values  mean  inflow/outflow. 


Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

Soya 

-0.6 

-0.7 

-0.9 

-1.1 

-1.0 

-1.1 

-1.4 

-1.2 

-0.8 

-0.7 

-0.5 

-0.5 

Tsugaru 

-1.5 

-1.4 

-1.5 

-1.5 

-1.3 

-1.3 

-1.5 

-1.5 

-1.2 

-1.5 

-1.4 

-1.6 

Tsushima 

2.1 

2.1 

2.4 

2.6 

2.3 

2.4 

2.9 

2.7 

2.0 

2.2 

1.9 

2.1 

d.  Optimal  determination  of  vertically  integrated 
velocity 

The  variational  problem  is  the  determination  of 
(Utj,  V,j,  A through  minimizing  the  cost  function  L, 


dL  i)L  i)L 

-  =  0,  -  =  0,  -  =  0. 

dUij  dVjj  d\  tj 

Substitution  of  (12)  into  (15)  leads  to 


(15) 


U„  =  Ujp  -  ^(A,  +  Ay+1  -  A;+li,  -  A;+1J+1)  (16) 
ViJ  =  V!f>  ~  2^(A"  +  K+1J  “  A-+1  ”  (17) 


+  u--'  ~  u'  w  ~  Wi) 

+  +  Vi-uj  -  Vu-1  -  K  y-  1  =  0.  (18) 

Substitution  of  (16)  and  (17)  into  (18)  leads  to  a  linear 
algebraic  equation  for  the  Lagrange  parameter, 

an^/-ij-i  5"  a2iA,j-i  T  c^Ay+jj-!  +  a12Xi_IJ  +  £122 Ay 

+  ct32Xi+ij  +  al3X ,-Lj+t  5"  a2iXjj+1  +  o3i\l+lj+1 


=  Sip  (19) 

where  i  =  2,  3,  •  •  • ,  Nz  —  1;  j  =  2,  3,  ■  ■  ■ ,  Ny  — 
1;  and 


a 


it 


a  22 


a2 1 

and 


Cly^  ^3i  U33 


1/J_  J_\ 

4\Ax2  +  Ay2)’ 


a  23 


a  12 


32  2  \  Ax2  Ay2 ) 


s.j  =  +  upix  -  un,j  -  un,j- 1) 


(20) 


+  ^(Vp  +  Vftj  ~  Vgli 


m^-i). 


(21) 


The  alternative  direction  implicit  method  (Press  et  al. 
1986)  is  used  to  obtain  the  value  of  the  Lagrange  pa¬ 
rameter  at  the  grid  point.  A,-,,  solving  (19).  Substituting 


the  given  value  of  A into  (16)  and  (17),  we  obtain  the 
optimal  estimation  of  U tj  and  Vij. 


e.  Determination  of  the  bottom  velocity 

Let  h  =  h(x,  y )  be  the  bottom  topography  and  (u  h, 
v_h)  be  the  bottom  velocity.  Applying  the  thermal  wind 
relation  (6)  to  the  two  levels  z  and  —h,  we  have 


(u,  v)  -  (i u ,  v)_h  =  J-  f  ly,  -y]  dz! .  (22) 

fpo  J-,,  \dy  dxJ 

Vertical  integration  of  (22)  from  the  bottom  (z  =  —h) 
to  the  surface  (z  =  0)  leads  to 

(«,  v)_h  =  \(U,  V ) 

li 


8 

fhPo 


(23) 


With  the  estimated  bottom  velocity  ( u ,  v)_h,  we  use  the 
thermal  wind  relation  (22)  to  obtain  the  absolute  ve¬ 
locity  from  the  density  field. 


f.  Volume  transport  streamfunction 


Due  to  the  continuity  (10),  the  volume  transport 
streamfunction  ("T)  can  be  defined  by 


U  =  - 


9y  ’ 


(24) 


and  satisfies  the  Poisson  equation 


/_cp  +  _  av  _  du 

ydx2  dy2)  dx  dy 


(25) 


5.  Absolute  geostrophic  velocity 

a.  Annual  mean 

Figure  5  shows  the  inverted  horizontal  velocity  vec¬ 
tors  at  depths  0,  50,  100,  150,  200,  and  300  m  respec¬ 
tively.  The  variational  P- vector  method  inverts  the  ve¬ 
locity  well  at  the  three  major  straits:  Tsushima/Korea, 
Tsugaru,  and  Soya  Straits.  However,  the  flow  in  Tatar 
Strait  is  not  well  resolved  due  to  the  poor  data  quality. 

We  take  the  velocity  field  at  depths  0,  50,  and  100 
m  (150,  200,  and  300  m)  to  represent  the  upper  (inter- 
mediate)-layer  circulation  features.  The  inverted  flow 
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Fig.  5.  Inverted  annual  mean  horizontal  velocity  vectors  at  different  depths:  (a)  0,  (b)  50  m,  (c) 
100  m,  (d)  150  m,  (e)  200  m,  and  (f)  300  m. 
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pattern  coincides  with  earlier  descriptions  of  the  JES 
circulation  (e.g.,  Uda  1934).  The  inverted  velocity  field 
agrees  well  with  early  results  using  the  P-vector  method 
(Chu  et  al.  2001a)  except  for  flow  in  the  straits,  which 
is  much  improved.  As  in  Chu  et  al.  (2001a),  the  NKCC 
along  the  Russian  coast  is  not  well  inverted  due  to  poor 
data  quality. 

1)  Upper  layer 

The  most  striking  feature  of  the  upper-layer  circu¬ 
lation  is  its  three-branch  structure.  North  of  35°N,  the 
TWC  bifurcates  into  an  eastern  portion  (first  branch, 
i.e.,  JNB)  and  a  western  portion.  The  strength  of  the 
TWC  at  both  portions  reduces  with  depth.  Flow  through 
the  western  portion  (i.e.,  EKWC)  closely  follows  the 
Korean  coast  until  it  separates  near  37°-38°N  into  two 
branches:  the  offshore  branch  (the  second  branch)  fol¬ 
lows  the  SPF  to  the  western  coast  of  Hokkaido  Island 
and  the  longshore  branch  (or  the  third  branch)  moves 
northward  (i.e.,  the  EKWC).  Such  a  three-branch  pattern 
was  first  identified  by  Suda  and  Hidaka  (1932)  and  Suda 
et  al.  (1932)  using  hydrographic  and  current  meter  data 
from  June  to  September  in  1929  and  in  1930,  respec¬ 
tively.  Since  then,  the  existence  of  three  branches  has 
been  believed  to  be  the  typical  TWC  flow  pattern.  We 
may  call  the  second  branch  the  Polar  Front  Current 
(PFC),  which  is  stronger  than  the  first  branch  (i.e.,  the 
JNB).  The  maximum  speed  of  the  PFC  is  found  around 
0.2  m  s  1  at  the  surface.  The  maximum  speed  of  the 
JNB  is  found  around  0.1  m  s-1,  occurring  near  the  Jap¬ 
anese  coast  (35.5°N,  135°E)  at  the  surface.  The  second 
feature  of  the  upper-layer  circulation  is  its  multieddy 
structure.  An  evident  cyclonic  gyre  is  identified  in  the 
Japan  Basin  (38°-44°N,  135°-140°E)  with  the  flow  as¬ 
sociating  with  the  PFC  as  its  southern  and  eastern  flanks 
and  the  flow  from  the  north  as  its  western  flank.  A 
mesoscale  anticyclonic  eddy  is  identified  in  the  UTB 
(36°-38°N,  130°-132°E)  with  the  PFC  as  its  northern 
and  western  flanks.  The  velocity  field  weakens  with 
depth. 

2)  Intermediate  layer 

The  most  striking  feature  of  the  intermediate-layer 
circulation  is  its  multieddy  structure.  The  PFC  weakens 
at  depth  150  m  and  disappears  at  depth  200  and  300  m. 
The  cyclonic  gyre  in  the  Japan  Basin  (38°-44°N,  135°- 
140°E)  is  evident  with  a  maximum  speed  around  0.1  m 
s_1.  The  UTB  anticyclonic  eddy  is  evident  at  depth  150- 
and  200-m  depth  and  weakens  drastically  at  300  m. 

b.  Seasonal  variability 

The  seasonal  variability  of  the  inverted  absolute  geo- 
strophic  velocity  (not  in  the  straits)  is  quite  similar  to 
an  early  study  using  the  ordinary  P-vector  method  (Chu 
et  al.  2001a). 


1)  Surface  circulation 

The  velocity  variability  reduces  with  depth  since  the 
seasonal  density  variability  weakens  with  depth.  The 
PFC  has  a  weak  seasonal  variation  in  flow  pattern  and 
a  strong  seasonal  variation  in  current  speed.  The  JNB 
is  too  weak  to  be  identified  in  the  winter  (Dec-Feb).  In 
March,  the  JNB  starts  to  occur  along  the  Japan  coast. 
It  strengthens  in  the  spring  and  reaches  the  maximum 
velocity  in  July. 

An  interesting  feature  in  the  surface  current  field  (Fig. 
6)  is  the  out-of-phase  variation  between  PFC  (along  the 
Polar  Front)  and  JNB  (along  the  west  coast  of  Japan). 
For  example,  the  PFC  (JNB)  strengthens  (weakens) 
from  July  to  September,  and  the  PFC  (JNB)  weakens 
(strengthens)  from  January  to  April. 

The  Japan  Basin  cyclonic  gyre  occurs  all  year  round, 
with  the  flow  associated  with  the  PFC  as  the  southern 
and  eastern  flanks  and  the  flow  from  the  north  as  the 
western  flank.  The  seasonal  variability  of  this  gyre  is 
largely  determined  by  the  seasonal  variability  of  the  PFC. 

2)  Out-of-phase  variation  between  PFC  and 

JNB 

The  seasonal  variation  of  the  JES  major  currents  is 
also  shown  in  the  meridional  cross-section  (135°E)  of 
the  monthly  mean  u  component  (Fig.  7).  Two  eastward 
flowing  currents,  PFC  (38°-40°N)  and  JNB  (36°-37°N), 
are  well  represented.  The  PFC  is  strong  from  September 
to  December  with  the  maximum  speed  of  0.2  m  s_1  and 
weak  from  February  to  June  with  a  minimum  speed  of 
0.1  m  s-1.  During  the  summer  monsoon  season  (Jun- 
Oct),  the  PFC  strengthens  (maximum  speed  increases 
from  0.12  to  0.18  m  s-1),  and  the  JNB  weakens  (max¬ 
imum  speed  decreases  from  0.12  to  0.04  m  s_1).  During 
the  winter  monsoon  season  (Dec-Apr),  the  PFC  weak¬ 
ens  (maximum  speed  decreases  from  0.18  to  0.12  m 
s  1 ),  and  the  JNB  strengthens  (maximum  speed  increas¬ 
es  from  0.04  to  0.12  m  s_1). 

3)  UTB  ANTICYCLONIC  EDDY 

The  seasonal  variation  of  the  UTB  anticyclonic  eddy  is 
shown  in  the  zonal  cross-section  (37°N)  of  monthly  mean 
v  component  (Fig.  8)  between  130°-132°E:  the  northward 
flow  in  the  west  and  the  southward  flow  in  the  east.  The 
inverted  UTB  eddy  has  a  core  from  the  surface  to  100-m 
with  a  maximum  swirl  velocity  of  0.06  m  s_1.  It  is  quite 
evident  during  the  whole  winter  season  (Dec-Feb)  and 
weakens  in  March.  It  cannot  be  identified  at  the  37°N  cross 
section  from  spring  to  early  summer  (Apr-Jun)  and  be¬ 
comes  evident  again  in  late  summer  (Aug). 

4)  Flow  through  the  Tsushima/Korea  Strait 

The  monthly  mean  surface  velocity  vector  field  in 
Tsushima/Korea  Strait  (Fig.  9)  shows  the  capability  of 
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Fig.  6.  Inverted  monthly  mean  surface  horizontal  velocity  vectors. 


the  variational  P-vector  method  to  invert  the  currents  in 
shallow  straits.  The  inverted  general  circulation  pattern 
agrees  with  the  observational  circulation  pattern  using 
acoustic  Doppler  current  profiler,  reported  by  Egawa  et 
al.  (1993):  The  main  axis  of  the  TWC  exists  in  the 


western  channel  at  34°N  for  all  months,  and  the  current 
in  the  eastern  channel  is  relatively  weak.  The  seasonal 
variation  of  the  current  magnitude  is  quite  small.  North¬ 
eastward  steady  flow  always  exists  in  the  middle  of  the 
strait  with  a  maximum  current  speed  of  0.1  m  s_1. 
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Fig.  7.  Meridional  cross-sections  of  the  inverted  monthly  mean  u  velocity  (cm  s~‘)  along  135°E. 
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Fig.  8.  Zonal  cross-sections  of  the  inverted  monthly  mean  v  velocity  (cm  s  ')  along  37°N. 
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Fig.  9.  Inverted  monthly  mean  surface  horizontal  velocity  vectors  in  Tsushima/Korea  Strait. 


6.  Volume  transport  through  the  Tsushima/Korea, 
Tsugaru,  and  Soya  Straits 

Warm  water  enters  the  JES  through  the  Tsushima/ 
Korea  Strait  with  the  TWC  from  the  East  China  Sea 
and  exits  the  JES  through  the  Tsugaru  and  Soya  Straits. 
There  is  no  evident  volume  transport  through  Tatar  Strait 
(Martin  and  Kawase  1998).  Using  the  inverted  absolute 
geostrophic  velocity  field,  we  calculate  the  monthly 


mean  volume  transport  through  the  four  straits.  The 
transport  through  Tatar  Strait  is  zero.  The  transports 
through  the  Tsushima/Korea,  Tsugaru,  and  Soya  Straits 
are  shown  in  Fig.  10.  An  evident  seasonal  variation  is 
found  in  Tsushima/Korea  Strait  with  a  minimum  value 
of  1.9  Sv  in  November  and  a  maximum  value  of  2.9 
Sv  in  July  (Table  2.)  The  annual  mean  volume  transport 
is  2.3  Sv,  and  the  range  (maximum  minus  minimum)  of 
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Fig.  10.  Monthly  mean  volume  transport  through  three  major 
straits:  Tsushima/Korea  (a:  solid),  Soya  (b:  dot-dashed),  and  Tsugaru 
(c:  dashed).  Here  positive  (negative)  values  indicate  inflow  (outflow). 

the  seasonal  variation  is  1.0  Sv.  This  range  (1.0  Sv)  is 
larger  than  Isobe’s  (1994)  estimation  (0.7  Sv)  using  the 
sea  level  difference  across  the  strait,  and  smaller  than 
Inoue  et  al.’s  (1985)  estimation  (1.61  Sv)  with  3.43  Sv 
in  summer  and  1.82  Sv  in  winter,  and  Yi’s  (1966)  es¬ 
timation  (1.9  Sv)  with  a  maximum  of  2.2  Sv  in  October, 
and  a  minimum  of  0.3  Sv  in  February. 

The  volume  transport  in  Soya  Strait  (outflow)  has  a 
similar  seasonal  variation  to  that  in  Tsushima/Korea 
Strait.  The  maximum  (minimum)  outflow  from  Soya 
Strait  is  1.4  Sv  (0.5  Sv)  in  July  (December).  The  annual 
range  is  0.9  Sv.  The  volume  transport  in  Tsugaru  Strait 
(outflow)  has  a  rather  weak  seasonal  variability.  The 
maximum  (minimum)  outflow  from  Tsugaru  Strait  is  1.6 
Sv  (1.2  Sv)  in  December  (September).  The  annual  range 
is  0.4  Sv.  The  ratio  of  the  outflow  through  Tsugaru 
versus  Soya  Strait  varies  from  a  maximum  value  of  3.2 
in  December  to  a  minimum  value  of  1.1  in  July. 

Table  2  shows  that  the  inverted  inflow  transport 
through  Tsushima/Korea  Strait  equals  the  total  outflow 
transport  through  Tsugaru  and  Soya  Straits  (a  mass  con¬ 
servation  pattern).  Such  a  feature  indicates  the  improve¬ 
ment  of  the  variational  P-vector  method  over  the  original 
P-vector  method,  and  might  be  practically  useful  in  de¬ 
termining  open  boundary  conditions  for  regional  nu¬ 
merical  models.  For  example,  some  JES  circulation  mod¬ 
els  use  constant  outflow  partitioning:  Chu  et  al.  (2001b) 
(or  Bang  et  al.  1996)  assumed  that  75%  (80%)  of  the 
total  inflow  transport  flows  out  of  the  JES  through  Tsu¬ 
garu  Strait,  and  25%  (20%)  through  Soya  Strait.  Such 
an  uncertainty  in  the  open  boundary  conditions  can  be 
eliminated  when  the  variational  P-vector  method  is  used. 

7.  Conclusions 

The  goal  of  this  study  is  to  investigate  the  seasonal 
variabilities  of  the  Japan/East  Sea  circulation  using  the 


U.S.  Navy’s  public  domain  Generalized  Digital  Environ¬ 
mental  Model  climatological  hydrographic  data  and  a 

recently  developed  variational  P-vector  method.  The  ma¬ 
jor  results  from  this  research  are  summarized  as  follows. 

1)  The  variational  P-vector  method  is  a  combined  box/ 
/3-spiral  model.  It  has  shown  the  capability  to  invert 
the  Japan/East  Sea  circulation  reasonably  well,  and 
in  particular  at  the  three  major  straits:  Tsushima/ 
Korea,  Tsugaru,  and  Soya  Straits.  The  inverted  ab¬ 
solute  geostrophic  velocity  field  coincides  with  ear¬ 
lier  observational  depiction  of  the  Japan/East  Sea 
circulation. 

2)  We  estimate  the  monthly  mean  transport  through 
three  major  straits:  Tsushima/Korea,  Tsugaru,  and 
Soya  Straits.  The  volume  transports  through  the 
Tsushima/Korea  and  Soya  Straits  have  evident  sea¬ 
sonal  variabilities.  The  annual  mean  volume  trans¬ 
port  through  the  Tsushima/Korea  Strait  is  2.3  Sv  with 
the  maximum  (minimum)  inflow  of  2.9  (1.9)  Sv  in 
July  (November).  The  annual  mean  volume  transport 
through  the  Soya  Strait  is  0.9  Sv  with  the  maximum 
(minimum)  outflow  of  1.4  (0.5)  Sv  in  July  (Decem¬ 
ber).  The  volume  transport  through  the  Tsugaru  Strait 
(outflow)  has  a  rather  weak  seasonal  variability.  The 
maximum  (minimum)  outflow  from  the  Tsugaru 
Strait  is  1.6  (1.2)  Sv  in  December  (September). 

3)  The  major  feature  of  the  inverted  annual  mean  cir¬ 
culation  is  the  three-branch  structure  in  the  upper 
layer  (surface  to  100  m)  and  the  multieddy  structure 
in  the  intermediate  layer  (100  to  300  m).  In  the  upper 
layer,  the  Tsushima  Warm  Current  bifurcates  north 
of  35°N  into  the  eastern  portion  [Japan  nearshore 
branch  i.e.,  the  first  branch]  and  a  western  portion 
[the  East  Korean  Warm  Current].  The  East  Korean 
Warm  Current  moves  northward  and  separates  near 
37°-38°N  into  two  branches:  the  offshore  branch  (the 
second  branch)  follows  the  Subpolar  Front  to  the 
western  coast  of  Hokkaido  Island,  and  the  longshore 
branch  (or  the  third  branch),  moves  northward.  This 
three-branch  pattern  weakens  with  depth.  In  the  in¬ 
termediate  layer,  an  evident  cyclonic  gyre  is  iden¬ 
tified  in  the  Japan  Basin  (JB;  38°-44°N,  135°-140°E) 
with  the  flow  associating  with  the  second  branch  as 
its  southern  and  eastern  flanks  and  the  flow  from  the 
north  as  its  western  flank.  The  maximum  swirl  speed 
of  this  gyre  is  around  0.1  m  s-1.  A  mesoscale  an- 
ticyclonic  eddy  is  identified  in  the  Ulleng/Tsushima 
Basin  (36°-38°N,  130°-132°E)  with  a  maximum 
swirl  speed  of  0.06  m  s~k  This  mesoscale  anticy- 
clonic  eddy  weakens  with  depth. 

4)  An  out-of-phase  variability  is  found  between  two 
major  eastward  currents:  the  Japan  nearshore  branch, 
and  the  Tsushima  Warm  Current  along  the  subpolar 
Front.  During  the  prevailing  summer  monsoon  sea¬ 
son  (Jun-Oct),  the  Tsushima  Warm  Current  along 
the  subpolar  Front  strengthens  with  the  maximum 
speed  increasing  from  0.12  to  0.18  m  s_1,  and  the 
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Japan  Nearshore  Branch  weakens  with  the  maximum 
speed  decreasing  from  0.12  to  0.04  m  s~k  During 
the  prevailing  winter  monsoon  season  (Dec-Apr), 
the  Tsushima  Warm  Current  along  the  Subpolar 
Front  weakens  with  the  maximum  speed  decreasing 
from  0.18  to  0.12  m  s  1 ,  and  the  Japan  nearshore 
branch  strengthens  with  the  maximum  speed  increas¬ 
ing  from  0.04  to  0.12  m  s-1. 

5)  The  mesoscale  anticyclonic  eddy  in  the  Ulleng/ 
Tsushima  Basin  has  seasonal  variability.  The  eddy 
is  quite  evident  during  the  whole  winter  season 
(Dec-Feb)  and  weakens  in  March.  It  cannot  be  iden¬ 
tified  at  the  37°N  cross-section  from  spring  to  early 
summer  (Apr-Jun)  and  appears  again  in  late  summer 
(Aug). 

6)  Currently,  the  JES  numerical  models  are  generally 
initialized  with  either  observed  or  climatological 
temperature  and  salinity,  as  well  as  zero  velocity. 
This  leads  to  an  imbalance  between  the  density  and 
velocity  fields  during  the  spinup  period.  The  Navy’s 
GDEM  T,  S  and  the  inverted  velocity  fields  provide 
balanced  initial  fields  for  Japan/East  Sea  numerical 
modeling  and  simulation. 
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